Interaction of a graphene sheet with a ferromagnetic metal plate 
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Nanoscale surface forces such as Casimir and the van der Waals forces can have a significant 
influence on fabrication, handling and assembly processes as well as the performance of micro and 
nano devices. In this paper, the investigation and the calculation of the Casimir force between 
a graphene sheet and a ferromagnetic metal substrate in a vacuum are presented. The reflection 
coefficients of graphene are graphene-conductivity dependent, and the conductivity of graphene is 
described by the Kubo formalism. There is an effect of magnetic properties of the metal on the 
Casimir interaction. The magnetic effect plays a significant role at low temperatures or high value of 
chemical potential. The numerical results also demonstrate that the thickness of a metal slab has a 
minor influence on the Casimir force. The investigation and findings about the Casimir force in this 
study would lead to useful information and effective solutions for design and manufacturing of micro 
and nano devices, especially in the areas of micro and nano machining, fabrication, manipulation, 
assembly and metrology. 

PACS numbers: 78.20.Ls, 31.30.jh, 78.67.Wj, 12.20.Fv 



I. INTRODUCTION 

Micro and nano technologies are potential economic 
engines and are considered as enabling technologies with 
exceptional economic capabilities. There have been in- 
creasingly more products in which the integrated micro 
and nano-electromechanical systems (MEMS/NEMS) 
are crucial in today's market because they allow for 
improved functionality, lower costs and higher quality. 
With the decrease in the size of MEMS and NEMS, addi- 
tional nanoscale surface forces, such as the Casimir force 
and the van der Waals force, should be considered [TJ [S], 
especially in the areas of micro and nano machining, fab- 
rication, manipulation, assembly and metrology. There- 
fore, a fundamental understanding of the Casimir force 
has recently been under intense discussion in the research 
and technological development (RTD) communities. 

The Casimir force can cause the small elements in a 
device to stick together. At small scales, the nanoscale 
surface forces may overcome elastic restoring actions in 
the device and lead to the plates' sticking during the 
fabrication process [Tj. It has been recently shown that 
Casimir forces may hamper the functioning of MEMS 
and NEMS devices by providing a pull-in instability [4]. 
The Casimir force between two objects can induce adhe- 
sion and stiction leading to failures in devices [H 0H2]- 
Particularly, when the size of a material is under the 
threshold length, the influence of the loading interaction 
on design and manufacturing of micro and nanodevices 
becomes much more prominent. 

The Casimir interaction originates from the quan- 
tum electromagnetic fluctuations between two objects 
[3] [6]. It provides a fundamental understanding 
about nanoscience and has a significant influence on 
the high performance of devices. Over six decades, the 



Lifshitz theory has been employed to investigate this 
force in metal-metal [7J [5] , semiconductor-semiconductor 
[S], metal-superconductor [TOJ [TT], metamaterials- 
metamaterials O [13] and graphene-graphene systems 
[Tmi7j . Results of some configurations have been ex- 
perimentally examined and have good agreement with 
theoretical calculations [HI [TH] . 

The interaction between a non-magnetic metal with a 
magneto-dielectric material has also recently been dis- 
cussed [T^I21j . These studies illustrate that the mag- 
netic properties can significantly influence in the Casimir 
pressure. Theoretical and experimental studies also show 
that it is possible to obtain the repulsive Casimir force 
in such systems [H [501 and the results could pro- 
vide possible solutions to handle the stiction and adhe- 
sion problems found in MEMS/NEMS devices. In terms 
of theory, the calculations depend substantially on the 
model that describes the dielectric function of metals. 

Graphene, an atomically thin layered material with 
novel properties, has gained a great deal of attention 
since its discovery [25] . Researchers from a wide range 
of scientific fields have spent considerable effort in ex- 
ploring its nature and its potential for practical ap- 
plications [HH27J. Recently obtained results demon- 
strate that graphene can be a promising candidate for 
next-generation electronic devices [28]. The functions of 
graphene have been applied to diverse areas from biology 
to material science. 

In this paper, we calculate the Casimir force between 
a graphene layer and a ferromagnetic (FM) substrate 
shown in FigJT] The FM material is taken to be Fe de- 
scribed via a dielectric and magnetic response properties. 
We utilize the Lifshitz theory to examine the role of the 
separation, temperature, and thickness of the substrate. 
We also investigate the influence of the magnetic response 
of the FM and show that in most cases its contribution is 
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= 2irkBTl/h is the Matsubara frequencies. The re- 
sponse properties of the metal are characterized by the 
dielectric function £i(i£z) and the permeability function 
fJ>i(i£l) which are frequency dependent along the imagi- 
nary axis {oji — We use the Drude model E\d and 
the plasma model E\p to describe the dielectric function 
of metals OH [2T] 
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FIG. 1: (Color online) Schematics of a two-dimensional layer 
(graphene) and a substrate (FM metal) with thickness D. 



relatively small compared to the one from the dielectric 
response. 

This paper is organized as follows: Sec. II presents 
the detailed expressions and calculations of the thermal 
Casimir interactions. The results and discussions are 
mentioned in Sec. III. Finally, Sec. IV presents the con- 
clusions. 
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where cu p is the plasma frequency and the damping pa- 
rameter is 7 P . At room temperature, the I = term of the 
Casimir-Lifshitz force in the metallic system is dominant. 
This reason allows us to consider the contribution of the 
static magnetic permeability /x(0) > 1 on the Casimir 
force. For higher orders of /, one should put = 1 

[21] . Both models lead to r^ M (0) = 1. However, we can 
obtain two different expressions of r^ E (0) when applying 
these two models (Hj 



II. CASIMIR INTERACTION BETWEEN 
GRAPHENE AND A METAL SUBSTRATE 

Here we calculate the Casimir force between the 
graphene sheet and a substrate in vacuum at a separa- 
tion a. The force per unit area at temperature T is given 

as 0115112211321 
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the reflection coefficients corresponding to the transverse 
magnetic (TM) and transverse electric (TE) field modes. 

The reflection coefficients for the substrate are given 
as gl|T2l|22| 
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It is important to note that fcj_ is the wave vector com- 
ponent perpendicular to the plate, c is the speed of light, 
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here D (0,k±) and P (0,k±) are the TE reflec- 
tion coefficients corresponding to the Drude and plasma 
model, respectively. 

For graphene, the reflection coefficients are found to 
be [1611221 
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where a = cr(i£,) is the 2D conductivity described via the 
the Kubo formalism [I|>] 
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The first term in Eq.([7| corresponds to the intra- 
band contribution, while the second term corresponds 
to the interband contribution. At low temperatures, the 
graphene conductivity a(i£i) approaches to the universal 
value cr = e 2 /4h. 

Using Eq.lJfjJ and Eq.Q the TM reflection coefficient 
of graphene is given by 
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FIG. 2: (Color online) The normalized graphene conductivity 
a(i£)/ao vs frequency at given temperatures. 
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where Xt — hc/kgT is the thermal wavelength in most 

( 2) 

materials. The explicit expression of r^E^O, fcj_) shows 

(2) 

that it is susceptible to temperature. r>p^(0,k±) is 
nonzero only because of the intraband component of the 
graphene conductivity. 

The first Matsubara frequencies £i are approximately 
0.027, 0.162, and 0.325 eV corresponding to the tem- 
perature 50, 300 and 600 K, respectively. As seen in 
Fig. [2j for / > 2, it is possible to substitute cr for <j(i£i) 
in Eq.Q to calculate the reflection coefficients and the 
higher order terms in the Casimir-Lifshitz force formula. 



III. RESULTS AND DISCUSSIONS 

In this section, we consider the Casimir interactions 
between graphene and a Fe metal plate. The parameters 
of Fe are uj p = 4.09 eV, j p = 0.018 eV, and /i(0) = 



10 4 [7]. As discussed in previous part 
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As a result, the I = term is the dominating term. Now, 
to calculate the contribution of the TE mode F^(a,T) 
with I = 0, it is necessary to choose a good model for Fe. 
Both the plasma [TH] and Drude model [3T] , however, has 
been widely used to compute the Casimir interaction and 
fit with experimental data. Determining an accuracy of 
two models as compared to measurement data has been 
a controversial issue. In the following work, we utilize 
the plasma model to calculate the interactions. The ex- 
pression of F^(a,T) is expressed by 
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Note that F^(a,T) is a unique term affected by the 
magnetic property of the ferromagnetic substrate. In 
Fig^J the ratio F^ ] E /F^ ] M is less than 0.75 % for dis- 
tances a in the range from to 1 \im at given tempera- 
tures. It indicates that the effect of /i(0) on the Casimir 
force can be neglectable. For this reason, Fe can be 
treated as a regular non-magnetic metal lx{iQ = 1- 

For a finite-thickness iron slab D, the reflection co- 
efficients of the metal material are modified as follows 
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The expression above is identical to the term with I = 
in the Casimir force between two metals. The sum of the 
I > 1 is much smaller than that of the metallic system 
due to the presence of the low graphene conductivity [TB] . 



To study the influence of thickness on the Casimir 
force, we rewrite Eq.|l]) in the following way 

F(a,T) = Ffy(a,T) + F®(a,T) +F^(a,T), (13) 

here (a, T) is the sum of all I > 1 terms and 
is thickness-dependent. Because Tj/ M (Q, fcj_) = 1, so 
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Wp M (0, k±) = 1. From this, the expression of F^^{a, T) 
is in Eq.( |10[ ) and is independent of D. The component 
force F^JJa, T) is weaker than that for the thick plate 
and one can consider F^ E {a,T) ps 0. Therefore, only 
F^(a.T) is sensitive to a variation of subtrate thick- 
ness. 
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FIG. 4: (Color online) The Casimir pressures corresponding 
to the approximate and full expression at different tempera- 
tures. Here F (a) = -ir 2 he /240a 4 . 

Figure [3j shows that the thickness of the Fe slab does 
not influence much the Casimir forces. At the same tem- 
perature, the separation between the two curves corre- 
sponding to the finite plate and the semi-infinite plate 
is small. In bulk materials, a reduction of the thick- 
ness gives rise to a small decrease of the Casimir force. 
The thermal wavelength of graphene is At/200 ~ 38 nm 
at room temperature. At distances beyond the ther- 
mal wavelength of a system, the contribution of the 
sum of I > 1 terms can be ignorable. In the investi- 
gated range, the influence of F^(a,T) is not significant. 
This finding suggest that since graphene is a two dimen- 
sional monolayer material, it significantly interacts only 
with the nearest interface layers of the bottom substrate. 
Graphenc's effect on layers separated by a sufficient num- 
ber of layers is negligable. The Fe slab with D = 20 nm 
has more layers than that which are meaningful. Another 
interesting special feature is that at large distance limit 
(a > 0.8fim) the contribution of F^\a,T) to the total 
Casimir force is minor. It can be explained that The 
/ > 2 terms provide weak force components due to the 
presence of the small universal conductivity in the reflec- 
tion coefficients. It is important to note that a universal 
graphene conductivity is a main reason for less than 2 % 
difference as one compares the dispersion force between 
two graphene sheets with Fo(a) [16]. It demonstrates 
that the approximate formula for the Casimir interaction 
equivalent to the I = term of the Casimir-Lifshitz ex- 
pression in this regime is F app (a,T) = — k B T(^(3)/8ira 3 . 



In Fig[4j we plot the normalized Casimir force with the 
full expression (see Eq.|l])) and the approximate expres- 
sion. The I = term is dominant and can replace the 
full expression of the Casimir interaction when a > 0.8 
/im for 300 K. 

The critical T c of Fe is 1043 K. At T > T c , the mag- 
netic properties of the ferromagnetic material nearly van- 
ish. Spins in the Fe slab are rearranged so that = 1 
at all frequencies. Nevertheless, the dielectric function 
is not affected by the directions of spins. Therefore, the 
plasma frequency and damping parameter are unchanged 
by temperature increases at T > T c . If we consider the 
Casimir interaction between Fe and another metal, the 
phase transition leads to a significant change in the mu- 
tual interaction HI] . The fact that graphene is quite 
transparent with the presence of the magnetic properties, 
the dispersion force may not be modified. 

A notable point is that the formula in Eq.Q was used 
for pure graphene with the chemical potential [ic = 0. 
If fic 7^ 0, the graphene conductivity has to be rewritten 
in the form of [33] 
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The first two terms are known as the intraband con- 
ductivity, another is the interband component. Eq.(14) 



shows that the intraband contribution of cr(/ics*£) still 
induces a non-zero value of r TE (0,k±). Obviously, 
the graphene conductivity is strongly susceptible to the 
chemical potential. The chemical potential of a graphene 
sheet can be controlled by using an applied electric field 
E d 021 [351 



ireoh 2 Vp 



E(f(E)-f(E + 2^ c ))dE, (15) 



here f(E) is the Fermi distribution function, Vf = c/300 
is the Fermi velocity. 

Other ways to modulate fic are an applied magnetic 
field [35] and chemical doping. Using a magnetic field also 
varys the expression of graphene conductivity and cre- 
ates Landau energy levels. In the limit of extremely low 
chemical potential, r!^ E (0,k±) is represented as Eq.jj)]). 
For He 3> k B T, the expression of the TE reflection coef- 
ficient with I = is written by 
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In Fig[5] the Casimir interactions are calculated at low 
temperatures (< 50 K) and fic — 1 eU with and with- 
out the magnetic properties. Unlike the case of a pristine 
graphene sheet (fic = e ^0> the influence of permeabil- 
ity of Fe on the dispersion force is noticeable at 10 K. 
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FIG. 5: (Color online) The Casimir forces at the case with 
and without taking into account the magnetic properties at 
He = 1 eV. 

The properties, however, have much less effect on the 
Casimir force at temperatures greater than 50 K. The 
majority of F(a,T) comes from the interband conduc- 
tivity of graphene. At 10 K, the contribution of the 
intraband on the Casimir force can be considerable to 
that of the interband conductivity. However, as tem- 
perature increases the interband contribution increases 
substantially, and causing a decrease in the contribution 
from intraband. For this reason, the magnetic properties 
nearly disappear at T > 50 K. 

To investigate further the Casimir force dependence 
on chemical potential, we consider the fluctuation inter- 
actions at a = 100 nm and various temperatures versus 
He shown in Fig{6] F(a,T)/ F changes from 0.018 to 
0.034 at T = 10 K and from 0.045 to 0.064 at T = 300 
K when \ic varys in the regime below 1.2 eV. In addi- 
tion, at higher temperatures, the curves of F(a,T)/ Fq is 
of the linear form. 

It is remarkable that Fe is not a superconducting mate- 
rial, so iron has no temperature phase transition as tem- 
perature decreases to 10 K. Alternatively, the universal 
value of conductivity <7o also demonstrates that even if 
temperature is K, graphene still does not change phase 
and properties. 

IV. CONCLUSIONS 

When design and manufacturing of micro and nano de- 
vices and components, it is important and necessary to 
understand the loading and effects of related forces on 
the systems and its components, especially the nanoscale 
surface forces such as Casimir and the van der Waals. In 
this study, a comprehensive investigation and discussion 
of the Casimir interaction between graphene and FM ma- 
terials is presented. We discussed about the Casimir force 




nj(eV) 

FIG. 6: (Color online) The Casimir pressures as a function of 
chemical potential. 

which is a function of the graphene conductivity and the 
permeability constant of Fe. It was shown that the con- 
ductivity of graphene heavily depends on the chemical 
potential which again has a large impact on the disper- 
sion force and causes a difference between non-magnetic 
and magnetic calculations at the ultra-low temperature. 
In other cases, there is no effect of the magnetic proper- 
ties on the Casimir force. One can apply a bias electric 
field to a graphene sheet to tailor the chemical poten- 
tial of graphene in order to control the magnitude of 
the Casimir force. According to our numerical results, 
the force between a graphene layer and the semi-infinite 
metal is smaller than the interaction between two metals 
in [21] because the conductivity of metals is higher than 
that of graphene. The numerical results also demon- 
strated that the thickness of a metal slab had a mi- 
nor influence on the Casimir force. The investigation 
and above findings about the Casimir force in this study 
would lead to useful information and effective solutions 
for design and manufacturing of micro and nano devices 
where graphene and FM metal materials are utilised, es- 
pecially in the areas of micro and nano machining, fabri- 
cation, manipulation, assembly and metrology. 
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